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Abstract 

We study the angle-resolved spontaneous emission of near-infrared light sources in 3D photonic 
crystals over a wavelength range from 1200 to 1550 nm. To this end PbSe quantum dots are used as 
light sources inside titania inverse opal photonic crystals. Strong deviations from the Lambertian 
emission profile are observed. An attenuation of 60 % is observed in the angle dependent radiant 
flux emitted from the samples due to photonic stop bands. At angles that correspond to the 
edges of the stop band the emitted flux is increased by up to 34 %. This increase is explained 
by the redistribution of Bragg-diffracted light over the available escape angles. The results are 
quantitatively explained by an expanded escape-function model. This model is based on diffusion 
theory and adapted to photonic crystals using band structure calculations. Our results are the 
first angular redistributions and escape functions measured at near-infrared, including telecom, 
wavelengths. In addition, this is the first time for this model to be applied to describe emission from 
samples that are optically thick for the excitation light and relatively thin for the photoluminesence 
light. 



I. INTRODUCTION 



There is much interest in nanophotonic control over the spontaneous emission of light 1 ' 2,3 . 
Therefore, emission properties of light sources like atoms, molecules, and quantum dots em- 
bedded in myriad nanostructures are studied intensively. The spontaneous emission rate of 
these light sources is described by Fermi's golden rule. The golden rule shows that the spon- 
taneous emission rate is proportional to a local property called the Local Density of Optical 
States (LDOS) 4 . Indeed, recent experiments have demonstrated that it is possible to control 
the LDOS, and thereby the spontaneous emission rate, using photonic crystals 1 ' 4 ' 5 ' 6,7 ' 8 ' 9 ' 10 . 

A second important aspect of photonic crystals is their ability to strongly influence the 
directional emission spectrum from embedded light sources. These spectral changes are 
caused by optical Bragg diffraction 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 ' 23 ' 24 . Although there are many 
qualitative studies on light transport in photonic crystals, there exist only a few quantitative 
papers. One of the central concerns regarding light emission from photonic crystals is the 
angular distribution, or escape function of the emitted light over exit angles. Fundamen- 
tal insight in this escape probability distribution is needed to explain the photonic crystal 
properties 25 ' 26 , and the crystals' influence on the radiative emission of the light source. 
Surprisingly, the escape function is often disregarded, which may easily lead to misinterpre- 
tations of experimental results. For example, changes in the emission spectrum, measured 
from low-efficiency light sources inside photonic crystals may mistakenly be attributed to 
changes in the LDOS, which cause an enhancement or inhibition of the sources' spontaneous 
emission power 52 . Hence, detailed knowledge about the escape function and a comparison 
to the emission spectra from non-photonic crystals are needed to find the cause of changes 
in the measured emission spectrum 17,27 . 

Most experimental work on light emission in photonic crystals and the corresponding 
escape function is limited to the visible range. For experiments on very strongly interacting 
photonic crystals made of semiconductors such as silicon 28 ' 29 , near-infrared emitting light 
sources, in particular versatile quantum dots such as PbSe and PbS (A > 1100 nm) are 
important to avoid absorption' 53 . The ability to control spontaneous emission and light 
propagation at these frequencies is of keen interest for applications, as spontaneous emission 
is a limiting factor in numerous devices 1 . Emission control therefore implies a substantial 
increase in efficiency in applications such as light beaming for the telecom industry, reflector 
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design, and solar cell light control 30 . Furthermore, guiding of light can be used to enhance 
the collection efficiency in single photon on demand applications. In this paper, emitted light 
escaping from 3D photonic crystals is studied at near-infrared wavelengths in the range of 
1200 to 1550 nm. To this end colloidal PbSe quantum dots are used as light sources inside 
strongly-photonic titania inverse-opal photonic crystals with a range of lattice parameters. 

II. EXPERIMENTAL DETAILS 

The samples consist of commercially available PbSe quantum dots (Evident Technologies) 
embedded inside titania inverse opal photonic crystals. See References 31,32 for a detailed 
description of the fabrication and characterization of these crystals. X-ray diffraction and 
scanning electron microscope studies show that the direction perpendicular to the samples' 
surface corresponds to the (111) direction of the crystal 32 ' 33 . 

The quantum dots emit at near infrared wavelengths A ~ 1300 nm, and we estimate the 
ratio between the lights' transport mean free path / and a typical sample thickness L to be 
l/L w 3 at these wavelengths 25 . Braggs' law was used to derive the required radii r of the 
polystyrene spheres that yields the crystals' lowest order stop band (which corresponds to 
the L-gap in band-structure calculations) at that emission frequency, taking into account 
a typical 30 % shrinking of the crystal during fabrication 32 ' 34 ' 35 . As a result, polystyrene 
sphere radii in the range of r = 340 to 620 nm are needed. Additionally, r < 200 nm spheres 
are used for reference measurements. For these radii, the lattice parameter is too small 
to cause Bragg diffraction at the quantum dot emission wavelength, and therefore these 
reference crystals are referred to as non-photonic. 

Not all fabricated structures are suited for emission experiments. Parts of the crystals 
may be covered with bulk titania, or the air-sphere ordering may be poor. Visually, samples 
that show clear Bragg diffraction are selected for further study. Subsequently, bright-field 
microscope images are made to chart each sample. Finally, reflectivity measurements are 
used to reveal the areas with clear stop bands that are suitable for subsequent emission 
experiments. 

The useful samples are infiltrated with colloidal PbSe quantum dots; similar to Refer- 
ences 5 ' 15 ' 17 ' 20 . Each sample is placed in a 5 ml glass vial. Next, 50 [A of 0.92 /iM PbSe 
quantum dots in hexane suspension is added. After one day the suspension is removed and 
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Figure 1: (a) Photograph of the sealed chamber used in the optical experiments. A transparent 
cylinder is sealed off by two plates. The bottom plate is mounted on four stages which are fixed 
to a rotation stage. The top plate contains a valve used to purge the chamber with nitrogen, (b) 
Front-view, gray-scale image of a titania inverse opal photonic crystal glued to a needle-pin (on 
the right). The white spot on the sample is a A = 532 nm laser spot that is used to locally excite 
PbSe quantum dots inside this structure. Due to a combination of diffusion, scattering and camera 
exposure time, the size of the laser spot appears much larger than the size of the actual focus 
(radius « 7/xm). 

the samples are rinsed three times with hexane, each time for 20 s. Subsequently, the sample 
is left to dry for one day. The infiltration leads to a distribution of quantum dots on the 
titania surfaces inside the inverse opals. From the concentration of dots in the suspension 
and the volume of the air spheres with typical radii of 385 nm, the distance between adjacent 
quantum dots is estimated to be 120 nm. This is sufficiently large to avoid reabsorption and 
energy-transfer processes between the dots 36 . In the results section below it is deduced that 
the quantum dots are really inside the photonic crystal. 

To prevent photo-oxidation of the quantum dots, the sample preparation and handling 
was carried out in a high quality nitrogen-purged glove box (MBraun, LabStar). For optical 
measurements the samples were glued to the side of a needle pin, mounted in a sealed 
chamber, and kept under a 1.6 bar nitrogen atmosphere. Figure l.(a) shows the home- 
build chamber that can be mounted on a rotation stage via a threaded ring at the bottom. 
Between this ring and the chamber, stages are visible that can be used to tilt and translate 
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Figure 2: Schematic representation of the experiment. Light from a A = 532 nm pump laser is 
coupled into a single mode fiber and used to excite quantum dots in a sample, through objective 
01. Near-infrared emission from the dots is collected with objective 02 and sent to a spectrometer. 
A longpass filter (F) blocks the laser light. An achromat lens (L) focuses the collected light onto 
the entrance slit of the spectrometer. An additional dichroic mirror (DM) is used to guide the 
visible light to a camera for alignment purposes. The sample is rotated together with the pump 
beam such that angle Op is fixed and the emission collection angle ©e changes. Consequently, 
the excitation conditions are the same for every angle Op. The arrows indicate the propagation 
direction of the light. 

the chamber. Figure l.(b) contains an image of a sample glued to a needle pin. As only the 
thin side of the crystal is in contact with the needle, both the frontside and backside of the 
sample remain accessible for the optical experiments, discussed below. 

Figure 2 describes the basic concept of the experiment. The sample is illuminated by 
a slightly focused laser beam of wavelength A = 532 nm. The excitation spot radius is 
estimated to be about 7 fim. Light emitted from the sample is collected by a 0.05 NA 
objective and sent to a spectrometer. An 850 nm longpass filter blocks the laser light in 
the detection path. The sample and excitation objective are rotated together, such that 
only the detection angle Ge changes, as in Reference 20 . After changing this detection angle 
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a CCD camera is used to roughly reposition the sample such that the illuminated spot is 
in focus with the spectrometer. Subsequent XYZ fine-tuning of the sample alignment was 
done by optimization to the emission signal measured with the spectrometer. The sample 
surface projected onto the entrance slit of the spectrometer changes with ©e- The measured 
intensities are corrected for this 54 , as discussed in Reference 20 . This correction also takes 
into account a phenomenological description of the excitation power distribution around the 
focal spot of the excitation beam. The current configuration allows the emission detection 
over an unmatched range of both positive and negative angles. 

Emission from the titania backbone was experimentally determined to be negligible at the 
emission wavelength of the quantum dots, whereas backbone emission did exert influence on 
the measured signals in earlier experiments at visible wavelengths 5 . Reasonably, the titania 
emission is expected from defect luminescence in the middle of the titania band gap, i.e., in 
the visible, and not around wavelengths A = 1300 nm. 

To reduce systematic errors, the emission detection angle 6e is varied in the following 
way. First, Oe is increased from 0° to maximum positive angles in steps of 20°. Subsequently, 
the angle is decreased in steps of 20° to maximum negative angles, which results in a data 
set with steps A0e = 10°. For some measurements the procedure is repeated to obtain a 
data set with steps A6e = 5°. The experiment ends where it started, at Oe = 0°. At the 
start and end of the experiment the excitation laser is blocked after which the background 
signal is determined. As a result, the reproducibility of the measurement can be judged 
for all frequencies, and in all experiments. For the typical integration times of 30 to 100 s, 
error bars of 1 % of the experimental value can be used to estimate the uncertainty in the 
measurement. However, this uncertainty does not account for systematic errors caused by 
spectral changes that may occur for long excitation times, or alignment issues. When the 
emission spectrum changes with time, saw-tooth-like shapes appear if the emission intensity 
is plotted versus the external angle 9e- Therefore, the additional uncertainty in the results 
can be evaluated after a measurement. 
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Figure 3: Normal-incidence reflectivity spectra measured from four titania inverse opals. These 
samples were grown from templates with different sphere radii r sp h e re) namely: (a) 655 nm, (b) 
515 nm, (c) 403 nm, and (d) 180 nm (left ordinate). These radii correspond with lattice parameters 
a of 1351 nm (a), 1114 nm (b), 764 nm (c), and 367 nm (d). Stop bands due to the L-gap appear 
as peaks in the reflectivity spectrum, (e) Shows the emission spectrum of PbSe quantum dots in a 
hexane suspension (right ordinate). The reflectivity is not defined in a narrow band around 15,800 
cm -1 because of interference by a reference laser beam inside the spectrometer. 



III. RESULTS AND DISCUSSION 
A. Angle resolved emission 

Reflectivity measurements on four samples with different lattice parameters are shown 
in Figure 3. The reflectivity setup was described in References 37 ' 38 . Only measurements 
that show clear Bragg diffraction are shown. The emission spectrum of PbSe quantum 
dots is included in the same graph. The stop bands are overlapping or to the red of the 
quantum dot emission spectrum, and far to the blue of the emission spectrum for reference 
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measurements. Therefore, the fabricated samples cover the desired wavelength range. The 
quality of the selected samples is comparable to samples used in earlier experiments in the 
visible range 5 ' 17 ' 20 . The areas that show the most powerful Bragg diffraction are used in 
subsequent emission experiments. 

To verify that the quantum dots are really inside the photonic crystal, a non-photonic 
sample was tested. Figure 4. (a) shows the result of the angle resolved emission experiment. 
The background subtracted data are shown as open squares. Subsequently each data point is 
corrected for the angle dependent collection efficiency (black dots). The signal is normalized 
to its maximum value. The measured intensity profile is the same for positive and negative 
angles, which is an indication of the diffusive nature of the emitted light. In Figure 4.(b) the 
data are compared to models for diffuse volume emitters and isotropic surface emitters. The 
former model is discussed in Section III B below, whereas the latter model yields an angle- 
independent, constant radiant. The data clearly agree with the volume emitter model. 
This means that the measured quantum dots are really inside the photonic crystal. The 
contribution from emitters at the samples' surface can be neglected, as the intensity clearly 
tends to zero at 6e = 90°. 

The quantum dot emission is determined to be from the bulk of the crystal. Furthermore, 
the distribution of the light emitted from a non-photonic crystal shows the expected diffuse 
emitter profile. These two prerequisites being satisfied, we can start to the use photonic crys- 
tals and continue the emission experiments. Figure 5 shows extensive data sets measured on 
a photonic sample. Figure 5. (a) Shows the angle and wavelength dependent emission power. 
A stop band is observed where the emission power is reduced. Increase of detection angle 6e 
results in a shift of the stop band towards shorter wavelengths. To a first approximation this 
behavior can be explained by Braggs' law (black line) 34 . In Figures 5.(b,c) cross sections of 
(a) measured at constant wavelengths or angle show the angular and wavelength dependence 
of the emitted power. The photonic sample behaves very distinctly from a Lambertian emit- 
ter. The stop band causes a decrease of the measured emission intensity over an angular 
range of approximately 40°. In (c) the frequency of the emission intensity peak changes 
with angle. The peak is shifted to the blue if the stop band overlaps the red side of the 
spectrum, i.e., at small angles. At larger angles the stop band has moved to the blue part 
of the spectrum and the emission peak is shifted somewhat to the red (not shown). For the 
largest angle the peak shifts back to the frequency of the dots' actual emission maximum. 
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Figure 4: Angle-resolved emission intensity of PbSe quantum dots in a non- photonic titania inverse- 
opal. The lattice parameter a = 637 nm and the measurement is done at A = 1400 nm. (a) Open 
squares are the measured, background subtracted, intensities. Filled dots show the measured data 
corrected for the detected area. The intensity tends to zero as the angle approaches e = 90°. 
(b) Normalized emission intensity versus detection angle. Filled dots are the same as in (a). The 
dashed line is the expectation for a diffuse volume emitter. Solid line shows the expected behavior 
for an isotropic surface emitter. Clearly, our measurements are in correspondence with the diffusion 
model. 
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Figure 5: Measured, background subtracted and surface corrected PbSe quantum dot emission 
spectra from a photonic titania inverse opal with lattice parameter a = 1114 nm. (a) Emission 
intensity as a function of wavelength and measured angle Oe for a photonic sample. A clear 
suppression of emission intensity is observed. As an estimate for the wavelength of this attenuation 
the Bragg condition is plotted (black curve). The white areas at wavelengths A = 1200 and 1550 nm 
mark the boundaries of the quantum dot emission spectrum, (b) Cross sections of (a) at specific 
wavelengths shows clear signs of stop bands, (c) Cross sections of (a) at specific angles. The 
measured emission spectra are angle dependent due to Bragg attenuation of the light in the stop 
band of the crystal and diffuse scattering. In (b,c) stop band frequencies are indicated by the 
horizontal lines below the curves. In this experiment spectral changes with time were negligible. 



This experiment demonstrates the influence of the photonic crystals stop bands on the emis- 
sion spectrum of the quantum dots. This stop band attenuation of a part of the emission 
spectrum resembles experimental findings using dye molecules 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17,19,20,21 ' 23 , but 
was not observed with quantum dots before 20 . An expanded escape function model is used 
to quantitatively interpret the measurement. This is the topic of the next section. 

B. Escape function model 

The model used to interpret the measured emission intensities is based on diffuse light 
transmission through multiple scattering, opaque media 39 ' 40 ' 41 ' 42 ' 43,44 . A complete description 
is given in References 20,25,45 . Here, the most important relations are summarized to explain 
the essence of this model. 

Consider light emission from a sample. The emitted intensity // e ) that escapes the 
sample at external angles between 0e = cos~ 1 (/i c ) and cos _1 (/i e + d/i e ) with respect to the 
surface normal of the sample is written as 

/i e )d/x e = I to t(u)P(uj; /i e )d/i e . (1) 
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Here, /tot(^) that exits the sample through the detection face is the spontaneously emitted 
intensity at frequency uj. P(u>; /x c ) is a normalized probability distribution function that 
describes the distribution of photons over the available escape angles. We call P(u; /x e ) 
the escape function and use it to compare the model to experimental data. To determine 
P(u; fi c ) experimentally Equation 1 is rewritten such that 

P(u; He) = T ( x ■ (2) 

The numerator is given by the measured intensity and the denominator J t ot(^) is deter- 
mined by summing the measured I(u;fj, e ) weighted by 2tt sin(©E)d6E to approximate the 
integration over 2n solid angle. 

For diffuse wave transmission through random media the escape function is analytically 
derived, and given by 43 ' 46 

3 

P(uj; fa) = -/i c [r c (w) + ■ [1 - R D (u; /i;)]. (3) 

Here, cos~ 1 (/i i ) is the angle inside the sample, which will be related to cos~ 1 ( / u c ) using 
Snell's law 35 . i? D (cj;/i;) is the angle and frequency dependent internal-reflection coefficient. 
Finally, t c {uj) = z c {uj)/1{uj) is the extrapolation length ratio, defined by the extrapolation 
length z e {uj) and the transport mean free path of the sample. z e (u) is the distance from 
the samples' surface, at which the diffuse intensity extrapolates to 0. t c (uj) can be expressed 
as a function of the angle-averaged reflectivity of the sample boundaries -Rd(w) 42 ' 43 ' 44 , 2. e., 



( \ 2 



1-Rd(u) 



(4) 



Furthermore, it was shown how to evaluate this average reflectivity Rv(<ju) from the angle 
and frequency dependent internal reflection coefficient Rd(w, /ii) 43 ' 44 ' 46 . This results in 

- 3C 2 (a;) + 2CiM 

/?D( " ) = 3C 2 {u,)-2CM + 2 ' Wlth (5) 

C n (u) = / Rb(uj; /ii) dfii. (6) 
Jo 

Hence, the escape function from Equation 3 can be evaluated if Ry>(uj; nd is known. 
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So far the derivation of the model was general and derived for random photonic media in 
which light transport is diffuse due to multiple scattering. This isotropic model for P(u; /x e ) 
successfully described experimental findings on random media when Rq{uj] was modeled 
using Fresnel's law in combination with an effective refractive index. However, this does 
not contain the stop bands observed in photonic crystals. Hence, a different approach is 
needed to find an expression for Rb(uj] and expand the escape-function model to describe 
photonic crystals, a procedure pioneered by our group 45 . 

The reflectivity from titania inverse opals is mainly determined by simultaneous Bragg 
diffraction from (111) and (200) planes 47 . The internal reflection coefficient is therefore 
modeled as the sum of two Gaussian reflection peaks 20 , i.e., 



with angle-dependent peak reflectivities Ri^i^i), widths Acji j2 (/ii), and center frequencies 
^i,2(A i i) 55 - These center frequencies cannot be described by simple Bragg diffraction by 111 
and 200 planes, since the two diffraction conditions cross. At this crossing a large avoided 
crossing occurs 47 . Therefore, band structure calculations are used to model u;i j2 (/ii) 45 . For 
Ri^ifJ-i) and Aui^ifJ-i) normal incidence reflectivity data are used. Furthermore, it is as- 
sumed that these four parameters are smooth functions of /i;. This expanded escape-function 
model can now be used to describe the results on photonic crystals. 

C. Comparison experiment with escape function model 

Figure 6 shows the escape function P(u; /i e ) over large frequency ranges, derived from 
experiments and calculated. Figure 6. (a) is derived from three experiments on samples with 
different lattice parameters. The results from the different experiments are separated by 
the horizontal, gray lines. The central part (0.7 < a/A < 0.92) shows a strong stop band 
in the escape probability due to the L-gap. As expected, the frequency of the stop band 
increases with external angle. An enhanced escape probability is observed at a/A ~ 0.84 
for ©e ~ 50°. The enhancement of the escape probability stems from a redistribution of 
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Figure 6: Escape probability P(u>; 0e)- ( a ) Experimental results from three data sets with different 
lattice parameters a, separated by horizontal lines. Top: a = 1351 nm, middle: a = 1103 nm, and 
bottom: a = 830 nm were determined from reflectivity measurements and assuming an average 
refractive index n e ff = 1.155. The low-frequency region shows the expected, Lambertian-like 
behavior. The central frequency region shows clear inhibition of the escape probability due to 
Bragg diffraction. Regions of enhanced escape probability are observed around a/A = 0.85 and 1.0 
for Oe ~ 50° and 0° respectively, (b) Result from the expanded escape-function model. For 
a/A < 1.2 the model is in agreement with the experimental findings. The stop band caused by 
Bragg diffraction from (111) lattice planes is observed, together with the regions of enhanced 
escape-probability. The higher order Bragg wave diffraction is more pronounced in the model. 
It is recognized by the inhibited escape probability region that starts at a/A = 1.3 for 0e = 
and continues to lower frequencies at larger angles. White curves show the lowest 6 bands along 
the LU direction, using parallel momentum (fey) conservation at the surface to transform internal 
propagation angles into external propagation angles. The horizontal arrows indicate the frequencies 
of the cross sections used in Figure 7 below. The common gray-scale range is limited for clarity. 
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the light that is Bragg diffracted inside the photonic crystal, as discussed below. At lower 
frequencies (a/A < 0.7) there is no sign of a stop band and the emission escape probability 
decreases continuously with E . At high frequencies (a/A > 0.92) another region is observed 
with enhanced escape probability at small Oe- The escape probability rapidly decreases at 
larger angles. Dark areas at the largest frequencies point out the presence of inhibited 
escape probabilities caused by Bragg-wave-coupling from the (200) lattice planes 45 . For 
applications this map can be used to find the lattice parameter a that yields the desired 
angular dependent escape probability at the wavelength of choice. 

Figure 6.(b) shows the calculation from the expanded escape-function model. The calcu- 
lated escape probability shows the same characteristics as the experimental results. There 
is a clear inhibition of escape probability starting at a/A ~ 0.8, for Oe = 0°. The central 
frequency of this inhibition shows the same dependence on 6e as the experimental results. 
Another line of inhibited escape probability runs from a/ A = 1.33 at Oe = 0° towards 
a/A = 0.95 at E = 50°. Multiple Bragg-wave-coupling causes the escape probability to be 
inhibited for angles larger than 40° and frequencies 0.88 < a/A < 1.03. This result agrees 
well with our experimental results. The regions with the measured, enhanced escape proba- 
bility where also predicted by the model. Clearly, the expanded escape-function adequately 
describes the experimental data measured from these strongly photonic crystals. The white 
curves show the lowest 6 LU bands from a band structure calculation used to describe exper- 
iments in the visible regime 45 . All band energies are multiplied by 1.125 such that the lowest 
bands overlap with the measured L-gap edges. The multiplication factor does strikingly not 
vary per sample. The difference between experiments and band structure calculations are 
subject to further study: more information on the titania distribution in these large lattice 
parameter crystals is needed to verify the model that describes the unit cell. Furthermore, 
dispersion causes the refractive index to be different in the near-infrared wavelength regime. 
Nevertheless, this first comparison already shows a very good overlap between the calculated 
band structure and the escape function. 

To compare the results from the model and the experiments in more detail, cross sec- 
tions were made at the frequencies indicated by the arrows in Figure 6.(b). In contrast to 
Figure 6, the escape function is now plotted versus fj, e = cos(6e). In this representation 
the escape function for isotropic emitters tends towards a line, simplifying comparison with 
experimental results. Figure 7. (a) shows the escape function versus external angle at three 
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Figure 7: Cross sections of escape function at different reduced frequencies a/A. Symbols denote 
the measurements and lines stem from the escape function models, (a) Results in photonic range. 
The solid lines result from the expanded escape-function model, and the dash-dotted line results 
from the regular escape-function model. Escape function versus cos(0e) shows inhibition and 
enhancement of the escape probability with respect to the isotropic distribution (dash-dotted line). 
Inside (outside) the stop bands the escape probability is inhibited (enhanced) up to 60 % (34 %) 
with respect to the isotropic distribution. The frequency of the stop band is marked by the 
horizontal lines below the curves. The stop band moves to higher frequencies at larger angles. For 
the a/A = 1.00 case, the upward and downward pointing arrows respectively indicate the enhanced 
and inhibited escape probability with respect to the isotropic model, (b) Results in non-photonic 
regime. Circles and squares are measured on the same sample with lattice parameter a = 220 nm. 
The diamonds are data collected from a sample with a = 830 nm. Only for a/A = 0.17 the 
experimental data deviate from the model (open squares), as a result of an emission intensity 
increase with time. 
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different reduced frequencies in the photonic regime (a/A G {0.77, 0.87, 1.0}). The stop band 
is clearly visible as the region where the escape probability is reduced with respect to the 
isotropic probability distribution. Outside the stop band the escape probability is higher 
than the isotropic model such that the escape probability is conserved. This probability 
enhancement is caused by a redistribution of the light inside the photonic crystal: for a cer- 
tain angular range the light Bragg diffracts back into the sample, where it becomes diffuse 
and escapes the sample at angles away from the Bragg condition. The enhancement is most 
profound close to the stop band. Further away from the stop band the data converge to 
the isotropic case. The expanded escape-function model agrees very well with the exper- 
imental data. The dash-dotted line shows the regular escape-function model for isotropic 
emitters. The experimental data appear even more explicit than the model predicts. Com- 
pared to isotropic emission, inhibition of the emitted power as high as 60 % is observed 
in the photonic crystal. Furthermore, the photonic crystal gave rise to emission power en- 
hancements as high as 34 %. For a/A = 0.77 and 0.87 the inhibition and enhancement is 
observed over a cone with a top angle of 2 • 6 E = 60°, which corresponds to a solid angle 
of 2tt(1 — cos(6e)) = 0.84 steradian. The pronounced inhibition at a/A = 0.87, between 
cos(Qe) = 0.97 to 0.71, corresponds to a solid angle of 1.6 steradian. 

Figure 7.(b) shows the escape function P{u;jj, e ) cross sections measured from two non- 
photonic samples. Symbols denote measured data and the line is calculated from the regular, 
isotropic model. There is no stop band. The data at a/A = 0.15 and 0.17 are measured 
from the same sample. The data at a/A = 0.15 and 0.58 match the regular escape-function 
model. The a/A = 0.17 data deviate from this model, i.e., especially at small angles 
(cos(Oe) ~ 1, open squares) the measured data are above the isotropic model. This de- 
viation is explained by an observed emission intensity increase with time, at the blue side of 
the emission spectrum (a/A = 0.17). This intensity increase results in a saw-tooth shaped 
measured angle-resolved emission spectrum, because of the applied order in the measured 
angles. At 6e = 0° the result of the very first and last measurement was averaged. Due to 
this intensity increase the normalization used to calculate I tot (uo) in Equation 2 becomes less 
accurate. This causes the measurement to deviate from the model. The modulation of the 
measured escape function with external angle is much smaller at the red side of the spec- 
trum (a/ A = 0.15). Although the blue side of the emission spectrum may have changed with 
time the rest of the spectrum has remained practically constant. This clearly illustrates the 
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effect of inhomogeneous broadening in the quantum dot ensemble, i.e., the measurements 
at various wavelengths probe different subsets of quantum dots. Clearly, the regular escape- 
function model fits the escape-function data from non-photonic titania inverse opals very 
well over a very large range of reduced frequencies. Hence, we are able to model the escape 
function both from photonic and non-photonic crystals, using our new expression for the 
internal- reflection coefficient Rd(ou; fa) (Equation 7). 

D. Discussion 

Photonic crystals can be divided into (1) optically thick samples where the thickness 
L is much larger than the transport mean free path I (7/L <C 1), and (2) optically thin 
samples where the thickness is about the mean free path of light, or less (Z/L ^ 1). Inside 
optically thick samples light that is Bragg diffracted back into the sample, scatters off 
defects and becomes diffuse. Subsequently, this light can exit the structure at different 
angles. This results in the redistribution of the emission intensity over escape angles. In 
emission measurements this redistribution causes an emission intensity increase at the edges 
of the stop band that is clearly not related to effects in the local density of optical states; see 
also 45 ' 48 . In optically thin samples most emitted light escapes the sample ballistically. There 
is only little scattering in these structures. Bragg diffraction, however, may still decrease 
the emission intensity considerably in specific directions 13 . 

The samples used in this study have l/L 3, which means that for the quantum dot 
emission the sample is in the intermediate regime between being optically thick (I <C L) and 
beeing optically thin (/ ^ L). As a consequence a significant part of the light emitted by the 
dots need not become diffuse before escaping the photonic crystal. Still our model which is 
based on diffusion theory describes the experiments very well. Although the diffuse nature 
of the emitted light may be questioned, it is believed that the models' requirement for a 
distributed light source is fulfilled as (1) the excitation of the quantum dots is diffuse because 
the sample is optically thick for the 532 nm excitation light, (2) the absorption length of 
the light is about two orders larger than the sample thickness, due to the low concentration 
of quantum dots, and (3) the emission of the quantum dots randomizes the wavevector k. 

The good correspondence between the model and the experimental findings demonstrate 
the applicability of the escape function model to a new regime in photoluminesence experi- 
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Figure 8: Schematic representation of different regimes in photoluminesence experiments. The 
regimes are defined by the transport mean fee path of the emitted light l em and the excitation light 
l exc , with respect to the sample thickness L. The regimes marked with A and B have been studied 
before. A new regime C is defined where the excitation light is diffuse and the emitted light is on 
the edge of being diffuse or ballistic. The escape function model, based on diffusion theory, was 
succesfully applied in this new regime. 



ments. This regime is illustrated in Figure 8 and denned by the transport mean free paths 
of the emitted light (l em ) an d the excitation light (l e xc), both with respect to the sample 
thickness L. In Figure 8 the regime marked by "A", where l em , l exc ^> L, corresponds to op- 
tically thin samples. Regime "B", with l em , l exc <C L, corresponds to optically thick samples. 
In these regimes the escape function has been used succesfully. In this paper a new regime 
"C" is added with l exc <C L and l em ~ L. As discussed above, the escape funtion model, 
based on diffusion theory and expanded to photonic crystals, was succesfully applied in this 
new regime. 
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Strong photoluminescence modifications may easily lead to misinterpretation of experi- 
mental findings. There exist various examples where an attenuation in a small frequency 
range of the photoluminescence emission spectrum is explained by a change in the sponta- 
neous emission properties of the emitter 11 ' 16 ' 49 ' 50 ' 51 . In contrast, the attenuation is caused by 
Bragg diffraction that causes a non-isotropic escape probability. This effect can be under- 
stood from Equation 1: the emitters properties and the local density of states are contained 
in I tot (u), whereas the experimentalist measures I(u;/j, e ). In photonic crystals these spec- 
tra cannot be linked without considering the escape function P(uj; /i c ). Without the escape 
function it cannot be recognized if light is redistributed over directions that do not show 
Bragg diffraction. Naturally, lifetime measurements could solve the matter, but were not 
addressed in the articles mentioned above. 

The escape function from titania inverse opals was also studied using externally injected 
light 45 , and, in the visible, using embedded CdSe quantum dot emitters 20 . A strong redis- 
tribution of the light over exit angles was observed, similar to our results. An interesting 
difference is that in the CdSe quantum dot emission spectra, no stop bands or spectral 
changes were observed, even though the relative linewidth AA/A = 0.05 was comparable 
to the relative stopband width ^ m 0.1. In contrast, we observe spectral stop bands, simi- 
lar to earlier work on dyeg 11 ' 12 ' 13 . 14 .^, 16,17,19,20,21,23^ j n PbSe emission spectra with a relative 

linewidth AA/A = 0.147. 

IV. CONCLUSIONS 

We have studied the angular distribution of near-infrared spontaneous emission (1200 < 
A < 1550) from PbSe quantum dots inside 3D titania inverse opals. We have observed angu- 
lar redistribution of the light over the exit angle, caused by a combination of light diffusion 
and Bragg diffraction in the photonic crystals. The escape function was extracted from the 
experimental data and explained with a diffusion model, expanded to photonic crystals. We 
found a very good agreement between the measurement and the model. This is the first 
time for the model to be applied to the regime of diffuse excitation in combination with 
optically rather thin samples for the emitted light. Furthermore, an interpretation of emis- 
sion enhancement at the sides of the stop bands is given in terms of angular redistribution 
of emission intensities. In the photonic crystal an escape inhibition of 60 % is observed as 
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well as an enhancement of 34 %, both with respect to the escape distribution from non- 
photonic samples. Typically the emission is enhanced or inhibited over solid angles larger 
than 0.84 steradian. The data presented form the first experimental evidence of angular 
redistribution of quantum dot emission from photonic crystals at near-infrared wavelengths, 
including telecom frequencies. 
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